Coherence properties of an atom laser 
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S '• Abstract 

^71^ ' We study the coherence properties of an atom laser, which operates by 

^f^ • extracting atoms from a gaseous Bose-Einstein condensate via a two-photon 

rn , Raman process, by analyzing a recent experiment [(Hagley et al., submitted 

vQ I to Phys. Rev. Lett. (1999)]. We obtain good agreement with the experimen- 

^^ ' tal data by solving the time-dependent Gross-Pitaevskii equation in three 
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dimensions both numerically and with a Thomas-Fermi model. The coher- 
ence length is strongly affected by the space-dependent phase developed by 
the condensate when the trapping potential is turned off. 
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One of the most exciting prospects resulting from the Bose- Einstein condensation (BEG) 
of alkah vapors [|l]-^ is the possibihty of producing an intense, coherent, and directed beam 
of matter waves, i.e., an atom laser. Indeed, prototype atom lasers have already been 
demonstrated [^,|^ . Potential atom laser applications include time-and-frequency standards, 
atom holography, and nanolithography. A critical element in the operation of an atom laser 
is the "output coupler" by which atoms are coherently extracted from the condensate 0. 
The design of this element is key to controlling the properties of the atom-laser beam [0 . 
At least two output-coupler mechanisms have been demonstrated. Condensate atoms have 
been extracted by rf magnetic pulses MM and by two-photon Raman transitions |TDI]. A 



quasi-continuous atom laser was demonstrated recently g] by using a rapid-fire sequence 
of laser pulses each of which caused condensate atoms to undergo a Raman transition that 
transferred momentum while simultaneously changing their internal state so that they were 
not trapped by the magnetic potential. Earlier theoretical studies of the properties of atom 
lasers [0,|rT|, made no comparisons of theory and experiment. This paper examines the 



coherence properties of atom laser wavepackets by analyzing a recent NIST experiment fl^ 
which probes such properties by measuring the decay of the interference contrast of two 
overlapping wavepackets outcoupled from a sodium atom condensate and separated by a 
variable delay time At. 

A parent condensate with wavefunction ^o(r, ti) is prepared at time ti. In one experiment 
the harmonic trapping potential was left on all the time, and we take ti = 0. In another 
experiment, the trapping potential was turned off and the condensate allowed to expand 
freely for up to ti = 5 ms. A 100 ns standing-wave laser pulse was applied at time ti with a 
wavelength A^ = 589 nm, detuned 600 MHz to the red of atomic resonance. This first laser 
pulse diffracts the condensate [|l^ to make two wavepackets ipf moving in the z direction 
with momenta ±2p, where p = hk = h/X^. We consider only the +2k wavepacket, since 
the problem is symmetric. At time ^2=^1 + At, the wavepacket evolves to 

V^+(r, t2) = 0i(r - Az, t2)e'''''e-"-^^' , (1) 
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where Er = -1^ and m is the atomic mass. The slowly varying envelope function (/)i is 
initially just a copy of the parent condensate wavefunction with norm lap <C 1: 0i(r,ti) = 
ail)o{Y,ti). In the experiment |ap ~ 0.02, and the momentum spread of V'o is very small 
compared to 2hk. The first wavepacket moves Az = vAt in time t2 — ti, where v = 2hk/m 
is the group velocity (60 yum/ms). 

A second standing wave laser pulse at time t2 = ti-|-At creates a second set of wavepackets 
ip2 ■ The combined number of atoms in the +2k wavepacket is {\ipf + '4'2\'^)r, where the 
brackets imply an integration over spatial coordinates. This fast +2k wavepacket soon 
clears the slowly expanding parent condensate and later can be imaged experimentally. The 
number of atoms in the +2k wavepacket is proportional to the following contrast function 
C(ti, At), defined so as to vary between and 1: 

C(t„At) = -( Mr-Az,t2) ^^ 
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occuring in Eq. (2) relates two points separated by distance Az and provides a measure 
of the spatial and temporal coherence of both the parent condensate and the outcoupled 
wavepackets. In the hypothetical case that the moving packets are plane waves (£i = ^ = 
1), then r(ti, At) = cos{4ERAt/h) varies between +1 when the wavepackets remain in phase 
(4:EjiAt/h = 2nn or At = nr, where r = ^^ = 10 /is for Na atoms) and —1 when they 
are out of phase (At = (n + |)t). The two packets constructively and destructive interfere 
in these two respective cases, giving a contrast function C(ti, At) which oscillates between 
1 and on a time scale of At = r/2. Actual condensate wavepackets of finite Thomas- 
Fermi radius ztf in the z-direction will physically separate after times on the order of 
tTF = '^ztf/v, after which F ^ and C(ti, At) -^ 1/2. Thus C(ti, At) oscillates rapidly 
between 1 and when At <^ Itf and Az <^ ztf, and approaches | when At > Itf and 
Az > Ztf- We will see that when ti is long enough that significant phase modulation has 
developed across the condensate due to the nonlinear mean field, then C(ti, At) drops to ^ 
in a time short compared to t^F- 

In the NIST experiment, the harmonic trap had frequencies ^ of 14 Hz, 28/a/2 Hz, 
and 28 Hz in the x, y, and z directions respectively, and a mean frequency of ^ = 28/\/2 
Hz. If the parent condensate has 1,500,000 atoms, ztf{x) = 22 /im and txF = 740 fis. 
The characteristic time for developing phase modulation (i.e., momentum spread) across 
the condensate is l/ui=8 ms. The experimental ti varied from to 5 ms, and At from 
to around 500 /is. Wavepacket images were taken about 6 ms after ti, long after the fast 
wavepacket has cleared the stationary, slowly expanding parent condensate. The number of 
atoms in the ±2k wavepackets could be counted using such images. 

In order to provide a reference normalization to remove errors incurred by shot-to-shot 
fluctuations in the total number of condensate atoms, since a new condensate had to be 
made to measure the contrast for each delay time, the NIST experiment actually utilized 
a second pair of standing wave pulses to produce a new set of ±2k wavepackets. The flrst 
pulse of the second pair was applied at time ta = ti + 3ms, after the fast wavepackets from 
the pulse pair at (ti,t2) have moved away from the parent condensate. The second pulse of 
the second pair was applied at t4 = ts + At + r/2, where r/2 = 5/is. Thus, when At <^ Itf, 
the contrast function C2(t3, At + r/2) for the second pulse pair is exactly out of phase with 
the contrast function Ci(ti, At) for the flrst pulse pair. The experimental images separately 
determine the number of atoms in the ±2k wavepackets from the (ti,t2) and the (t3,t4) 
pulse pairs. The necessity to provide a normalization of the number of atoms each packet 
from shot to shot is avoided in the experiment by reporting as the "signal" the function 

^^'^' '^' ^'^ - C7,(t„At) + C2(t3,At + r/2) • (') 

Just like C, the signal S oscillates rapidly between and 1 for At <^ txF and approaches | 
when At > trF- We can deflne a "coherence time" Ate to be the time for the envelope of 
the signal function S to decay halfway from its At = value of 1 to its long time limiting 
value of ^, that is, S'(ti, ts. Ate) = 0.75. A corresponding "coherence length" is Azc = vAtc- 



The time-dependent Gross-Pitaevskii (TDGP) equation describes the dynamics of 
if) (r,t), which includes the parent condensate plus the fast ±2/;; wavepackets: 

dih Ti^ 

th-^ = -^ V'V^(r, t) + (Krap(r, t) + Vlaser(r, t)) ^{t, t) 

+ f/oiV|V^(r,t)|V(r,t), (5) 

where Uq = A-kTl a/m^ a is the s-wave scattering length, N is the total number of condensate 
atoms and Vl;rap(r,i) is the trapping potential. The contrast functions Ci and C2 can be 
calculated from tlj{r,t). The interaction of condensate atoms with the four standing-wave 
laser pulses can be written as 

4 
I1aser(r, t) = Vl COS (2ki ■ r) J2 fit - tn, St) , (6) 

n=l 

where VLf(t,6t) is the single laser pulse envelope. The laser pulse duration 6t = 100 ns is 
short compared to At. The factors t„ are the times at which the four experimental pulses 
are applied. 

The modification to the condensate wavefunction caused by a short-time, low-intensity, 
standing-wave laser pulse can be best understood in momentum space. Before the first pulse, 
there is only a component centered at p = 0. We can use time-dependent perturbation theory 
to show that sidebands appear after the pulse at p = ±2h'k. For our experimental conditions, 
it can be shown that these sidebands have the same shape as the p = component and have 
amplitude proportional to the pulse area. At short time e immediately after the n-th pulse, 

^{r, tn + e)^ V-olr, t„) - ^VLA{6t) 

2h 
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The pulse is applied to the parent condensate at time t„, n = 1,2,3,4, where the parent 
condensate wavefunction has evolved from ti to t„. The laser pulse area, VLA{6t), is the 
area under the single laser pulse appearing in Eq. (||), and ViA/h -C 1 in the experiments. 
After normalization, the wavepacket in Eq. (^ serves as the initial condition for subsequent 
evolution. 

We have used two different methods to evolve tp{r, t). The first is a numerical propagation 
of the three-dimensional (3D) solution to the TDGP equation, using a fast-Fourier-transform 
method with the slowly- varying-envelope approximation to reduce the size of the spatial grid. 
The slowly varying envelope approximation is excellent here because the momentum spread 
of ijjo is very small compared to 2hk. We have verified that this methodology, which will 
be described in detail elsewhere, gives excellent agreement with full numerical solutions of 
the TDGP equation in one- and two- dimensions. This allows us to calculate exact contrast 
functions for a zero temperature condensate for any time sequence of trapping potential and 
laser pulses. 

We have also calculated the contrast functions using a second approximate method which 
we call the time-dependent Thomas-Fermi (TDTF) method. Let us first consider the case 
when the trapping potential is turned off at t = prior to the first pulse at ti. Once 
Vi;rap is removed, the parent condensate, ipo{r,t), evolves freely, develops phase modulation 



and expands somewhat. The 3D form of 7/^0(1", t) can be easily found since, for expanding 
condensates where the TF approximation is vahd, the solution of the TDGP is self-similar, 
i.e., it can be transformed to its original shape (before release) by suitable axis scalings. 
The time dependence of the scale parameters has been shown [jl3| to obey coupled nonlinear 
ordinary differential equations. Once the atoms in high momentum states clear the p = 
condensate, they evolve as free particles (if VlA/Ti ^ 1) and move with velocity ±2/;,k/m. 
In our 3D model, the full condensate wavefunction thus evolves after the passage of the first 
pulse as follows: 

V;(r,t) ^ ^o(r,t) - ;^l^L^(5t)e-^^^«(*"*^)/'^ 
in 

X e2*'^^'-V^o(r-v(t-ti),ti) 

+ e-2*''^"-^o(r + v(t-ti),ii)]- (8) 

Using Eqs. (|^ and (H) we can develop the condensate wavefunction for any number of pulses 
and delays. When the trap is left on, the only modification to the above analysis is that the 
p = condensate does not develop phase or expand. 

Fig. |I] compares the calculated results for S'(ti,t3, At) with the NIST data [|T2[ for the 



case where the trap was held on. In each panel the experimental signal is plotted against the 
delay At used for the first pair of pulses. The signal was measured for At in 1 /xs increments 
up to At = 50 /is after which the increment was 30 /is up to At = 530 /xs. Fig. [l|a shows 
excellent agreement with the short time data, which was normalized to unity at the first 
peak at 10 /xs. The TDGP and TDTF calculations also agree well, except that the phase 
of the latter slightly lags that of the former because of the small acceleration of the fast 
wavepackets by the effective potential provided by the parent condensate. The long-time 
evolution of the signal envelope agrees very well between the TDGP and TDTF calculations. 
The coherence time Ate predicted by the two models, around 275 /is, is slightly longer than 
the measured value of 225 ± 40 /is. When the trap is on, the decay of S'(ti, At) is simply 
due to the reduction of the time-dependent overlap of the moving outcoupled wavepackets. 
Consequently, the calculated and measured coherence lengths, Az^ = wAtc = 17 /xm and 14 
/im respectively, are an appreciable fraction of the size of the parent condensate, ztf = 22 
/im. This also implies that coherence extends across most of the size of the outcoupled 
wavepackets. This result is in line with a recent estimate of the coherence length of a static 



condensate using Bragg spectroscopy [14 



Fig. |] compares the experimental data for 5'(ti,t3, At) with the TDGP and TDTF cal- 
culations for two cases for which the trap was turned off at t = 0. In Fig. ^ (ti,t3) = 
(1.2ms,4.2ms), whereas in Fig. ^, (ti,t3) = (5ms, 8ms). The agreement between the two 
calculations, as well as the agreement between experiment and theory, is good for both 
cases. The coherence times and lengths are much smaller for these trap-off cases than for 
the trap-on case in Fig. |I[ For Fig. ^a the respective TDGP and TDTF At^ are 82 and 80 
/is as compared to 65 ± 10 /is for the experiment. For Fig. ^ the corresponding theoretical 
values of 38 /xs and 37 /xs compare with a measured value of 45 ± 10 /is. The respective 
coherence lengths for the (1.2ms,4.2ms) and (5ms, 8ms) cases are 5 /xm and 2 /xm, much 
smaller than ztf- Since Azc is substantially smaller than the condensate size, there must 
be another source of coherence loss than wavepacket separation. 



The extra source of coherence loss when the trap is off is due to the particle interac- 
tions that give rise to the nonlinear term in the GP equation. When the trap potential is 
removed, the parent condensate experiences the effective potential iVt/olV'oP) which causes 
phase modulation to develop across the condensate. This is due to the increased spread 
in the condensate momentum distribution as the atoms accelerate. For example, Fig. 1 of 



Ref. |I^ shows the spatial oscillations in Keipo and Imipo due to this phase modulation. The 
presence of these oscillations in iljQ{r,t) spoil the phase matching when packet 1 is trans- 
lated by Az during the interval At, and lead to a much faster loss of coherence between the 
packets than for the trap-on case. The longer ti is, the greater the coherence loss will be. 
Since the characteristic time scale to reach terminal momentum spread is 1/uj = 8 ms, much 
coherence loss is to be expected for the example in Fig. §3. 

In conclusion, outcoupled wavepacket coherence times and lengths predicted by a three- 
dimensional time-dependent Thomas Fermi model are in excellent agreement with those 
calculated from solution of the three-dimensional TDGP equation and also give good agree- 
ment with data from a recent experiment which measured coherence properties of Raman- 
outcoupled atom laser wavepackets. Since the outcoupled wavepackets are copies of the 
parent condensate, the experiment probes both the coherence of the parent condensate as 
well as that of the outcoupled wavepackets. Spatial and temporal coherence is maintained 
across the parent condensate while the trap is left on, but is rapidly lost when the trap is 
turned off due to phase modulation which develops across the condensate. 
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FIGURES 

FIG. 1. Comparison of calculated TDGP (solid line) and TDTF (dashed line) and experimental 
(points) signal functions S{ti = 0,1^ = 3ms, At) versus At for the case where the trap with 
1,500,000 atoms was held on during the laser-pulse firings, (a) Comparison during the first 50 ^s 
where the delay was stepped in increments of 1 fis. (b) Comparison of the TDGP signal with the 
measured signal envelope over the full delay range to 500 fis. The TDTF model gives essentially 
the same envelope. 

FIG. 2. Comparison of TDGP (solid line) and TDTF (dashed line) signal functions S{ti,t3, At) 
with the data (points) for the case where the trap potential was turned off at t = 0. (a) (ti,t3 = 
(1.2ms, 4.2ms) for 500,000 atoms in the trap, (b) (^1,^3) = (5ms, 8ms) for 2,500,000 atoms in the 
trap. 
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